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Abstract 

Extending the standard model with three right-handed neutrinos 
{Nk) and a second Higgs doublet (r?), odd under the discrete parity 
symmetry Z2, Majorana neutrino masses can be generated at 1-loop 
order. In the resulting model, the lightest stable particle, either a bo- 
son or a fermion, might be a dark matter candidate. Here we assume 
a specific mass spectrum (Mi <^ M2 < M3 < m^) and derive its con- 
sequences for dark matter and collider phenomenology. We show that 
(i) the lightest right-handed neutrino is a warm dark matter particle 
that can give a ~10% contribution to the dark matter density; (ii) 
several decay branching ratios of the charged scalar can be predicted 
from measured neutrino data. Especially interesting is that large lep- 
ton flavour violating rates in muon and tau final states are expected. 
Finally, we derive upper bounds on the right-handed neutrino Yukawa 
couplings from the current experimental limit on Br{fi — > 67). 



1 Introduction 



Solar [T], atmospheric [2] and reactor [3] neutrino experiments have demon- 
strated that neutrinos have mass and non-zero mixing angles among the dif- 
ferent generations. On the other hand observations of the cosmic microwave 
background, primordial abundances of light elements and large scale struc- 
ture formation have firmly established that most of the mass of the Universe 
consists of dark matter (DM) These experimental results are at present 
the most important evidences for physics beyond the standard model. 

There are several ways in which neutrino masses can be generated. Cer- 
tainly the best-known mechanism to generate small Majorana neutrino masses 
is the seesaw [S]. However, a large variety of models exist in which lepton 
number is broken near-or at- the electroweak scale. Examples are supersym- 
metric models with explicit or spontaneous breaking of R-parity [6], [7], models 
with Higgs triplets [8], pure radiative models at 1-loop [9] or at 2- loop [10] 
order and models in which neutrino masses are induced by leptoquark inter- 
actions [TT] . 

According to their free-streaming length DM particle candidates can be 
classified as either hot, warm or cold DM. Due to their large free-streaming 
length the mass and density of hot DM particles are strongly constrained [12)H 




Contrary, cold DM particles have a free-streaming length which is irrelevant 
for cosmological structure formation. Actually, cold DM is usually consid- 
ered the best choice to fit large scale structure data [13]. Warm DM (WDM) 
particles, for instance those that decouple very early from the thermal back- 
ground, have a smaller temperature than that of hot dark matter relics and 
thus a shorter free-streaming length. 

It has been argued in the literature [12] that WDM scenarios may be able 
to overcome the shortcomings of the standard cold DM scenario. Constraints 
on WDM particles have been quoted in reference [16] . If DM consists only of 
WDM, 

''^wDM ^ 1-2 keV whereas in mixed scenarios, in which the DM relic 
density receives contributions from cold and WDM as well, m^DM ^16 eVE|. 

The question of whether neutrino mass generation and DM are related has 
lead to a large number of models [T7]. In this paper we focus on a particular 
realization, namely the radiative seesaw model [18]. In this scheme three 
right-handed neutrinos Ni and a second Higgs doublet r] = {rj'^.rp), odd 

^Models where all dark matter is hot are ruled out completely by current cosmological 
data [13] 

^In this case WDM gives a contribution of ~10% to the total DM relic density [16] 
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under the discrete parity symmetry Z2, are added to the standard modeL 
As a result (a) the new Higgs doublet has a zero vacuum expectation value 
and there is no Dirac mass term. Thus, neutrinos remain massless at tree 
level; [b) the lightest particle in the spectrum, either a boson or a fermion, 
is stable and therefore, in principle, can be a dark matter candidate [20] . 

Here we study the implications for DM and possible collider signatures of 
this model. Our analysis is done in a particular scenario in which the Yukawa 
couplings of N3 are larger than those from and the right-handed neutrino 
spectrum is such that Mi <C M2 < M3. The right-handed neutrinos are 
assumed to be always lighter than the charged and neutral scalars. As it will 
be shown the lightest neutrino singlet can not be a cold DM candidate and 
instead behaves as WDM, contributing with less than 10% to the total DM 
relic density. In addition we will show that current experimental neutrino 
data enforces a number of constraints on the parameter space of the model. 
These constraints, in turn, can be used to predict the decay patterns of the 
charged scalar 77^. Therefore, the hypothesis that this model is responsible 
for the generation of neutrino masses (within our scenario) and that A'^i is a 
WDM particle can be tested in collider experiments. 

The rest of this paper is organized as follows: in section [2] we briefly 
describe the model, paying special attention to the neutrino mass generation 
mechanism. In section [3] we present simple and useful analytical results 
for neutrino masses and mixing angles. In section H] we discuss dark matter 
within the model and show that the lightest right-handed neutrino is a WDM 
relic. We then turn to the collider phenomenology of charged scalars in 
section [51 We show that different ratios of branching ratios of r/^ can be 
predicted from measured neutrino mixing angles. In section [6] we analyse the 
implications of the model for lepton flavour violating decays, in particular 
for fi 67. Finally in section [7| we present our conclusions. 

2 Neutrino mass generation 

The model we consider [TH] is a simple extension of the standard model, 
containing three SU(2)l x U{1)y fermionic singlets Ni and a second Higgs 
doublet rj. In addition, an exact Z2 discrete symmetry is assumed such that 
the new fields are odd under Z2 whereas the standard model fields are even. 
The Yukawa interactions induced by the new Higgs doublet are given by 

C = e,tKjNjPLLlr]' + h.c. (1) 
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Figure 1: Feynman diagram for Majorana neutrino masses. A = R,I labels 
the contributions from the neutral CP-even and CP-odd components of the 
Higgs doublet rj. 

Here, L are the left-handed lepton doublets, a, j are generation indices (Greek 
indices label lepton flavour e,fi,T) and eat is the completely antisymmetric 
tensor. Apart from these Yukawa interactions the quartic scalar term 

Ih<I>v)\ (2) 

where is the standard model Higgs doublet, is also relevant for neutrino 
mass generation. Since Z2 is assumed to be an exact symmetry of the model 
7] has zero vacuum expectation value. Thus, there is no mixing between the 
neutral CP-even (CP-odd) components of the Higgs doublets. The physical 
scalar bosons are, therefore. Re 0°, 77^, 77^ = Ker]^ and 77° = Imrj^. 

The setup of equation ([1]) and equation ([2]) generates Majorana neutrino 
masses through the diagram shown in figure [H The resulting neutrino mass 
matrix can be written as 

ca Mk hak hpk Bo{0,m\,Ml). (3) 

^ A=I,R 
k=1...3 

Here A = R,I, Mk are the right-handed neutrino masses, rriA are the 77^4 
masses, cr = +1 while c/ = —1 and Bo{0,m\, M^) is a Passarino-Veltman 
function [19]. The function Bq has a finite and an infinite part. Note that 
the infinite part cancel after summing over A and the resulting formula can 
be expressed as a difference of two Bq functions. The finite part of the 
Passarino-Veltman function B^ is given by 



As pointed out in reference [TS] if rj^ and rji are almost degenerate, i.e. 
m|, — mj = 2\^iP (t;2 = {2^/2Gf)~^) is assumed to be small compared to 
ml = (m|. + mj)/2, the neutrino mass matrix in ([3]) can be rewritten as 

Depending on the relative size between mo and M^ this formula can be 
simplified [18]. Here we will focus on the limiting case m^ ^ Ml. 



3 Analytical results 

Here we will consider a right-handed neutrino spectrum such that Mi <^ 
M2 < M3. In addition, as previously mentioned, we will also consider the 
limiting case ^ M|. In this case the neutrino mass matrix in eq. ([5]) 
becomes 

A 

®^ "^0 fc=1...3 

In general the neutrino mass matrix receives contributions from diagrams 
involving the three right-handed neutrinos. However, if A'^i is light enough, 
let us say, 0{Mi/M2) < 10~^ the contributions from A^i become negligible. 
In this limit Det[7Vli,] ~ and therefore only two neutrinos have non-zero 
masses. In this case simple analytical formulas involving neutrino mixing 
angles and Yukawa couplings can be derived. Note that in this limit only a 
hierarchical spectrum is possible. In what follows we will focus on the normal 
spectrum. Some comments on the inverted one will be given in section 15.21 
In the limit Det[A1,y] ^ the mass matrix structure is determined by 
the Yukawa couplings ha{2,3)- Therefore, it is useful to define two vectors in 
parameter space 

ha = {hi2, h22, h32) , 

hs = (/il3,/i23,/i33) ■ (7) 

In terms of these vectors the two non-zero neutrino masses can be written as 
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where ^/ is given by 



and 



rN 



(r-iv|h2 

^ M2 
" Ms" 



(9) 



(10) 



The ratio between the solar and the atmospheric scale is approximately given 

by 



R 



m 



1^2 



m 



(11) 



1^3 



Thus, from eq. ([HD and (Q, it can be noted that R is independent of ^/ and 
therefore independent of A5 and mo- 

The generation of the non-zero lightest neutrino mass can be under- 
stood from the misalignment angle between the parameter space vectors h2,3 
(cos 6' = h2 ■ ha/ 1 112! I ha I) which, from eq. ([HD, can be written as 



sin^^ 



1-R 
1 + R 



(12) 



where hr = |h2p/|h3p. Note that since hr as well as are positive quanti- 
ties a complete alignment between h2 and hs (sin^ = 0) is only possible if 
R = 0. However, this possibility is excluded as it implies m,y^ = 0. 

There is a minimum value of sin^ 9 consistent with the experimentally 
measured values of R. This value is determined by 



sin^^l 



{l + K r^f 



R 



1 + R 



7D 



1 + Rn^ir 



(13) 



and corresponds to the minimum misalignment between h2 and hs. Thus, 
in order to reproduce the correct solar and atmospheric mass scale ratio 
sin^ ^ > 0.47. Figure [2] shows the misalignment allowed region. 

Although not consistent with neutrino experimental data there is an in- 
teresting limit when the contribution from N2 to the neutrino mass matrix 
is small in comparison with those from A^^s. In this case the neutrino mass 
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Figure 2: Allowed range of the misalignment between the vectors h2 and h^. 

matrix becomes projective and therefore it can be diagonalized with only two 
rotations. The rotation angles can be written as 

tan6'23 = -T— , 

tan^i3 = -^=^==. (14) 

V "-23 + '^33 

As it will be shown in section [5] these results are good approximations in the 
case we are considering. 

4 Fermionic dark matter 

Before discussing possible collider signals of the charged scalar [f] we will study 
the implications of the assumed mass spectrum, Mi -C M2 < M3 < m^, 
on DM. In ref. [20], A^i was assumed to be a cold DM particle. Based on 
this assumption, two crucial observations, related with and the Yukawa 
couplings hai, were made: 

(i) The following relation has to be satisfied in order to obtain the observed 
■^We will denote the 77* mass by m,. 
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DM relic density, Qah^ ^0.12 [1]: 



Restricting the Yukawa couplings to the perturbative regime, i.e. the 
left-hand side of f fTSj) < 1, it was found that < 350 GeV. Fur- 
thermore, the constraint f|T5|) . being a lower bound for the Yukawa 
couplings hai, should be compared with the constraint derived from 
/i — > 67, which gives an upper bound for the Yukawa couplings (see 
section El eq. (142|) ). The apparent contradiction between these bounds 
was overcome by assuming a specific structure for the Yukawa couplings 
in reference [20] . 

(ii) The constraint Mi > 10 GeV for > 100 GeV must be satisfied in 
addition to the requirement that Mi < m^. 

If (i) and (ii) are combined, the hierarchical mass relation M1/M2 < 0(10"^) 
imply that M2,3 > which is not consistent with the analysis of neutrino 
masses discussed in the previous section. Moreover, this relation, in turn, re- 
quires another suppression mechanism for /i — 67 Q Therefore, the assumed 
mass spectrum. Mi <C M2 < M3 < m^, does not fit within the cold DM 
scenario of |20j . 

In what follows we will discuss whether A'^i can be a viable WDM candi- 
date. In this case there are important differences compared with the conven- 
tional sterile neutrino WDM scenario in which sterile neutrinos are produced 
by non-resonant active-sterile neutrino oscillations [221 123], namely: 

(a) The decay of Ni is forbidden by the Z2 symmetry [f]. Thus, the X- 
ray constraint [2^ [25l [26l [27] . derived from the absence of detection 
of X-ray photons from sterile neutrino radiative decays, can not be 
applied. This constraint, when applied to the conventional sterile neu- 
trino WDM scenario, yields an upper bound of ttIwdm ^ 4 keV [27] . 
This result combined with the Lyman-alpha forest data, which lead to 

suppression mechanism based on a low-energy flavor symmetry in the same type of 
models, with a radiative neutrino mass generation, was proposed in [3T|. 
possible origin of this symmetry was discussed in |17| . 
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a lower limit of m^DM ^ 10-14 keV, has ruled out the possibility |§ that 
all the DM consists of sterile neutrinos [SHIED] (see also [SI]). 

(b) In the conventional scenario the Yukawa couplings of the right-handed 
neutrino are tiny. Actually they can not be thermalized without mix- 
ing with the active neutrinos [22] and therefore can not be regarded as 
thermal relics. In contrast to the conventional case, the Yukawa cou- 
plings hak in the current model are not necessarily small (see section [6]). 
Thus, A^^i can be in thermal equilibrium at high temperatures. This im- 
plies that the constraints discussed in the literature on thermal WDM 
particles [121 [ISl [23 [30] can be applied in our case. Of course, the 
largest value of hak niust be consistent with the upper bound derived 
from — > 67 (see eq. (142|) ). 

Current cosmological data constraints [T2l [T6| [29| [30] imply that DM can 
consists of only Ni if the relativistic degrees of freedom at the decoupling 
temperature {g^{TD)) are larger than 10^, for Mi < 1 keV [16]. This is not 
satisfied in this model, the relativistic degrees of freedom can be at most 116. 
Therefore, A^^i can be regarded as WDM if there exists, in addition to A'^i, a 
dominant cold DM relic that gives a contribution of ~ 90% to the total DM 
relic density and if Mi < 16 eV [16] (this possibility, within the conventional 
WDM sterile neutrino scenario, has been throughout studied in |32j). 

From a more detailed analysis of this scenario we have found that the 
annihilation rate of A'^i at temperature T can be written as 

^ ^ ^ \ / ri a,f3 

Here we have assumed ^ T ^ Mi. The decoupling temperature can be 
estimated by equating the annihilation rate with the expansion rate, H = 
1.66 vCT^ rympi. From H{Td) ^ T{Td) we get 

^If the sterile neutrinos as WDM are generated in decays of some heavier particles, 
then the situation may change |28j . 
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where T^, is the decouphng temperature of the active neutrinos and g^{Ty) = 
10.75. For Td^2 GeV, for which g^iTo) = 77.5 0, eq. ([IT]) becomes 

/100GeV\^ , , , 

yii :^2.2xl0-^ 18 

\ nir, J 



which, as we can see from eq. ( I42l) . satisfies the constraint coming from 
H —>■ ej. Note that a stringent experimental upper hmit on Br{fi — ej) 
will imply a larger decoupling temperature. For example, a three orders of 
magnitude more stringent bound on Br{fi — > 67), as the one expected in 
near future experiments [33], will enforce To to be larger than ~ 140 GeV. 



5 Collider physics 

The Yukawa couplings that govern neutrino physics also determine the fermionic 
two-body decays of 77^ j and 77^. According to the Yukawa interactions in ([T]) 
possible decays of these states are: 

vh-'^aN, (19) 
V^^i^N,. (20) 

As will be discussed below A^2,3 follow decay chains that can lead to only 
missing energy. In that case the observation of the neutral Higgses 77^ j will be 
problematic. On the contrary, since charged scalar final states always contain 
-at least- a charged lepton their decays are easier to observe. Therefore, we 
will focus on charged Higgs decays. Apart from the Yukawa interactions the 
scalar doublet 77 has also gauge (and scalar) interactions which induce the 
decays rj^ ^ Tj\j W^, if kinematically possible. 

At LHC charged scalars can be produced either in association with a 
neutral scalar (single production) or in pairs [31]. In the former case the 
mechanism proceeds via qq annihilation mediated by a virtual W vector 
boson whereas in the latter case through s-channel exchange of a virtual 7 
and a Z: 

m ^ v^vh (21) 

qq 7/"^7/~ . (22) 
^We have assumed that at Td only iVi, among iV^'s, remains relativistic. 
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Charged scalar production in association with an 77° has been calculated in 
reference [3S]- According to this result the production cross section is larger 
than 100 fb for < 200 GeV. The pair production cross section, on the 
other hand, is further suppressed as it can not exceed 10 fb for charged 
scalar masses below 250 GeV [3lj. Contrary, at ILC the pair production 
cross section is larger than 10 fb for m^, < 350 GeV [2Sj- Thus, depending 
on the accumulated luminosity, LHC (ILC) should be able to explore up to 
masses of order ~ 200 - 250 GeV (400 GeV). 

5.1 Right-handed neutrinos: Decays, signals and iden- 
tification 

The correlations between charged scalar decays and neutrino mixing angles 
which will be discussed latter could be studied in collider experiments only 
if the decaying right-handed neutrino can be identified. Experimentally, in 
principle, this can be done. Let us discuss this in more detail: right-handed 
neutrinos, stemming from charged scalar decays, will produce, via an off-shell 
?7^, charged leptons through the decay chains 



In addition to these decay chains there are others which involve neutral 
scalars and lead to either dilepton + missing energy {i'^^^^a'^/s'Ni) or miss- 
ing energy {i^ai^/s^a'^is'Ni or I'a^/sNi) signals. 

The most important signatures for the identification procedure are (l23l) 
and (124|) due to their low backgrounds [37]. The right handed neutrino 
identification from the remaining decay chains might be rather hard as they 
involve additional missing energy. Thus, in general, they will diminish the 
relevant signals. Whether the decay branching ratios for the processes in 
f l23|) and fl2^ can dominate depend upon the different parameters (mainly 
Yukawa couplings and scalar masses), which we will now discuss in turn. The 
decay chains in eq. fl23l) dominates over the processes —>■ d-^tpi'a'J^fi'Ni 
and iVa — »• I'a^/s^a'^/s'Ni as long as 



(23) 
(24) 



J2 5r(iV3 ^ it^^u^'Up'N,) 



J2 ^KA^3 - f^^qC^^,>N,) > 



a',/3' 



a', 13' 



(25) 



a,/3 
a',/3' 
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The conditions on the parameter space of the model for which fl25p is ful- 
filled can be entirely determined from the three-body decay processes A'^j — > 
£^£jiVj and iVj — *• Vai-'pNj as the branching ratios in (125|) are given by 



a,l3 a,l3 
a', 13' a', 13' 

(26) 

Thus, from eq. (126|) and using the shorthand notation 

Br{N,^Nj) = Y,Br{N,-.i^i^N^), (27) 

a,/3 

Br,^,{N, Nj) = Y,Br{N,^u^upN,), (28) 

a.,13 

the constraints in fl25l) become 

5^(iV2 ^ iVi) >Br,UN2 ^ iVi) 
5^(iV3 _^ N2) X 5r(iV2 ^ iVi) >5nnv(iV3 ^ iV2) X 5ri„.(iV2 iVi) • (29) 

Similar conditions can be also obtained in the case of the decay chains in 
(!2l|) . namely 

Br{Ni^ Nj ) > 5ri,„ (N, ^ Nj) . (30) 

The partial decay width for the process Ni fafpNj summed over all 
possible final states for a fixed j, is given by 

Er(iv. - /./.iv.) = ' ' ' 38^3 ^m^m (31) 

where 

/(x) = 1 - Sx^ - 24x^ ln(x) + - (32) 

and S' = ?7if/ = £orS' = ?7r,/ if / = i^. This expression, in addition to the 
conditions (l29l) and (l30l) . lead to the constraint 



> . (33) 



Consequently, as long as the neutral scalars become heavier than the charged 
one the decay processes in 0231) and ([21]) become dominant. Note that this 
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result holds only if A^2,3 decay inside the detector. Whether this is indeed the 
case depends on the parameters that define eq. Since right-handed neu- 
trino masses M2,3 as well as the parameter space vectors I h2, 3] are bounded 
by neutrino physics, once the constraint fl33|) is imposed □ the only free pa- 
rameter is hi. Accordingly, the right-handed neutrino decay lengths are 
strongly determined by the value of |hi|. We calculate N2 and decay 
lengths by randomly varying the Yukawa couplings h^i for the benchmark 
point TJiRj = 140 GeV, = 150 GeV, M2 = 25 GeV and M3 = 45 GeV. 
After imposing neutrino physics contraints at la level [38] we get 

L2 C [0.08, 300] m, L3 C [10-^ 2] m , (34) 

which shows that A'3 always decay within the detector whereas N2 decays 
might occur outside. 

As can be seen from eq. (1311) the larger (smaller) |hi| the smaller (larger) 
L2. For the benchmark point we have considered, it has been found that 
in those regions of parameter space in which L2 is smaller than few meters 
Br{N^ —>■ N2) ~ (9(10~^) which implies that most decays will proceed 
through the decay chains in (12^ . On the contrary, when L2 is large N2 will 
behave, from the collider point of view, as A'^i and the only possible signals 
will be either dilepton + missing energy or missing energy. In this case 

Ni) > 0.1). 

In general, since from eq. fl3T]) we have 



according to our results the process N3 — s> i^^^^i ^^^^ ^e sizable {Br{N^ 



Br{N2 ^ iVi) = 4"^"^'' , (35) 

^R,I + ""^J? 

if -C mjii small values of |hi| will enhance the decays in (1251) . For the 
smallest value of |hi| for which A^2 still decays inside the detector (typically 
10"'^) we found that 

Br{N3 N2) X Br{N2 iVi) < 0.5 . (36) 

Hard leptons with missing energy (eqs. (125]) and (I24p ) are typical accel- 
erator signatures in conserving and non- conserving R-parity violating super- 
symmetric models [391 EZ] • Indeed, as pointed out in references [391 EZ] ? the 



^Scalar masses are also constrained from the requirement of scalar production at LHC 
or ILC (upper boimd) and from LEP data (lower bound). 
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discovery of supersymmetry could arise from such a signal. In the present 
case the possibility of having in addition displaced vertices might facilitate 
the reconstruction of N2 and iVs. Actually, since W and Z leptonic decay 
modes occur at the interaction point, this type of signals are practically back- 
ground free once the dipleton invariant mass distribution from the displaced 
vertex is above 10 GeV [37]. 

Regarding the identification procedure if decay according to fl23p the 
identification might be possible by counting the number of leptons emerging 
from a given vertex. In contrast to the decay chain (l23!l . if follows the 
processes in ( 12^ the number of leptons from decays will be the same, and 
the charged lepton counting "method" can not be used. In this case A^3 from 
A'"2 decays can be distinguished by looking to the kinematic endpoint of the 
lepton pair invariant mass distribution. This method have been extensively 
discussed in the MSSM context |10] and might be also applicable in this case. 
Note that the kinematic endpoint technique could be also applicable when 
follows the decay chain fl23l) . Thus, the right-handed neutrino identification 
procedure can be entirely based on this method. 

5.2 Collider signals related to neutrino physics 

The results presented below were obtained by numerically diagonalizing eq. 

for random parameters and checking for consistency with experimental 
neutrino constraints Different correlations among neutrino mixing an- 
gles and charged scalar decay branching ratios were found as expected from 
eq. (|T^ . The parameter mo, which essentially corresponds to mn or m/, was 
taken in the range 100 GeV < mg < 400 GeV[^ whereas the masses of A^3 and 
N2 between 40 GeV < M3 < 50 GeV and 20 GeV < M2 < 30 GeV 0. The 
Yukawa couplings were chosen such that |h2|/|h3| C [0.4,0.9]. In regions of 
parameter space in which N^^ and N2 are comparable -though lighter- to 
the correlations, discussed below, are less pronounced. However, the decays 
chains (see eqs. (l23l) and (!24l) ) will involve hard leptons from which the right- 
handed neutrinos can be readily identified. On the other hand, if A^3 and N2 
are much more lighter than rj^ the data points become strongly correlated. 
In this case, in contrast to the previous one, charged leptons emerging from 
the decay chains might be near the r -and possibly /x- threshold which will 

^The charged scalar mass was also taken in this range. 
^°Afi was taken below 16 eV as required by DM constraints 
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Figure 3: Ratio of decay branching ratios Br^± ^ / Br'^^^^ versus tan^ 6'23. The 
vertical strip indicates the current 3o" range for tan^ ^23 whereas the horizontal 
lines indicate the predicted range for this observable. 

render the right-handed neutrino identification problematic. 

Figure [3] shows that the ratio of decay branching ratios Br'^i''^ / Br^±^ 

{Br{ri^ Nk) = Br^J^'') is correlated with tan^^23- From the best fit 

point value (tan^6'23 = 1) Br^^'"^ ^ Br'^^^ is expected. Furthermore, the 3a 
range for the atmospheric mixing angle allows to predict this observable to 
lie within the interval [0.35,3.0], as indicated by the horizontal dashed lines 
in figure [31 

We have found that there exits an upper bound on the ratio of decay 
branching ratios 



strated by figure HI From this bound Br'^± ^ is expected to be around 1 order 
of magnitude smaller than Br^^^ + Br"^^^ . Which in turn implies, according 




(37) 

'13 and demon- 




to Br^r ^ SrT', that eA^s final states are further suppressed than jJiN^ 
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Figure 4: Ratio of decay branching ratios Br'^^±^ / {Br^± ^ + Br^^ ^) versus 
tan^6'i3. The vertical hne indicates the current 3a upper bound for tan^ 6'i3 
whereas the horizontal lines indicate the predicted upper bound for this ob- 
servable. 

and r final states. 

From eqs. ([8]) and (fTTj) we found a quantity, R^/R+, which is related to 
R. Rj: can be written as 



where i, j run over e, /i, r, tat corresponds to the right-handed neutrino mass 
ratio defined in eq. ( fTOl) and 7?./ is a common global factor that involves 
the same parameters that define Qf (see eq. (|9])) and decay branching ra- 
tios. Note that in the ratio -R_/i?+ this factor cancel. Numerical results 
are shown in figure O The spread in the plot is due to an ambiguity in the 
sign of the Yukawa couplings. From the current 3cr range for Ami2/Am23 
(vertical shaded strip in figure [5]) this quantity is predicted to lie in the range 




1 -4r^ 




, (38) 
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Figure 5: Ratio of decay branching ratios versus Ami2/Am23. The 

vertical shaded strip indicate the current 3cr range for Ami2/Am23 whereas 
the horizontal dashed lines shows the allowed region for 



(horizontal dashed lines) [1.4 x 10"^, 2.0 x 10"^]. 

As long as the constraints |h2|/|h3| < 1 and M2/M3 < 1 are satisfied 
the contributions of A^2 to the neutrino mass matrix are small in comparison 
with those from N^. Thus, the atmospheric and reactor angles approximate 
to eqs. ( 1141) . The results shown in figs. [3] and H] can be understood as a 
consequence of these constraints. Note that the sharpest correlations among 
the decay patterns of the charged scalar with neutrino mixing angles are 
obtained for the largest allowed (by neutrino experimental data) hierarchies 
between the parameter space vectors |h2| and jhal and the right-handed 
neutrino masses M2 and M3. 

In order to generate an inverted light neutrino mass spectrum {M.p)ii has 
to be large (of the same order of (A^;/)22,33,23)- Thus, large contributions from 
the loop involving N2 are necessary. These contributions spoil the leading 
projective nature of the neutrino mass matrix and therefore the atmospheric 
and reactor angles are no longer determined by eqs. f[T^ . Accordingly, the 
correlations among collider observables and neutrino mixing angles we have 
discussed will not hold in this case. However, in principle, these results can 
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Figure 6: Charged scalar fermionic decay branching ratio as a function of 
the average Yukawa couphng \h\ for the case in which W gauge boson final 
states are kinematically open. The charged scalar mass has been fixed to 
300 GeV. The shaded region is excluded by the experimental upper bound 
on Br{^ e'-f). 

be used to discriminate between the normal and inverted mass hierarchies as 
follows: If M3 > M2 and ^„ 5r^-^2 / Er^-^^ < i are experimentally 
established but none of the observables given in figs. [3], H] are found to be in 
the range predicted by neutrino physics the normal mass spectrum will be 
excluded. 

final states are also possible depending on whether the mass difference 
AM = — mjij is larger than M^j/ . Once kinematically open, this decay 
channel will dominate over the fermionic final states. However, as illustrated 
in figure [6], even in that case the fermionic decay branching ratios can be 
as large as ~ 10~^. Albeit possibly problematic to be measured at LHC 
might be measurable at ILC. As indicated in figure [6] (shaded region) larger 
values of these branching ratios are excluded by the current upper bound on 
i?r(yU — > 67) (see next section). 

^^This relation is derived from the constraint |h2|/|h3| < 1 
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6 Flavour violating charged lepton decays 



In this section we will derive upper bounds on the Yukawa couplings, hak, 
and briefly discuss possible low energy lepton flavour violating signals of the 
model. The set of Yukawa interactions induced by the right-handed neutrinos 
and the SU{2) doublet rj are responsible for lepton flavour violating radiative 
decays of the type la Ipl- Here we will concentrate on /i — > 67. The 
bounds derived from r — > 67 and r — /i7 decays are much more weaker than 
those from — > 67 and thus we will not consider them. 

In the limit mp -C rria the partial decay width of — //37, induced by 

and A^^, can be written as PI 



967r2 



^ IZlLi hakhpk\ 



nxt 



(39) 



From the above expression the decay branching ratio for /i — *• 67 can be 
written as 



Br {11 — s> 67) 



r(/i ^ 67) 



a 



247rG2, 



(40) 



and the current upper bound on this process yields the upper bound 



5^ Kkh2k 
k=l 



< 4.1 X lo- 



rn^ 



100 GeVJ 



(41) 



The largest value for these Yukawa couplings is derived from the largest 
charged scalar mass, m„ = 400 GeV, in this case 



fc=i 



Kk^^k 



< 6.5 X 10' 



(42) 



For smaller charged scalar masses the bound becomes more stringent. Note 
that this constraint can be satisfled by either hik <^ h2k or h2k <^ hik- An 
exception being the case of non-hierarchical Yukawa couplings. Under this 
assumption, and for = 300 GeV, the upper bound h < 1.9 x 10^^ can be 
placed. This constraint corresponds to the shaded region shown in flgure El 



""^^Note that we are considering the case <C mf^ and therefore the foop function 
reduces to a factor of 1/6. 
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Figure 7: Br{fi 67) as a function of the charged scalar mass under the 
assumption of non-hierarchical Yukawa couplings. The shaded region is ex- 
cluded by the current experimental upper bound. 

Since we do not have a theory for the Yukawas, an absolute value for 
Br{fi —>■ 67) can not be predicted. However, assuming non-hierarchical cou- 
plings this branching ratio is found to be larger than 10~^^ for h > 10~^ as 
shown in figure [71 Note that this result is a consequence of the assumption 
0{hik) ~ 0{h2k) and not a general feature of the model. 

7 Summary 

Assuming the mass spectrum Mi <^ M2 < M3 < we have studied some 
phenomenological aspects of the radiative seesaw model [18]. In particu- 
lar, we showed that current experimental neutrino data can be used to con- 
straint the parameter space of the model. Thus, some fermionic decays of 
the charged scalar 77^ are correlated with neutrino mixing angles which in 
turn allows to predict several ratios of decay branching ratios. Especially 
interesting is that if the 77^ is produced at colliders similar number of events 
with r and /i final states are expected, whereas events with e are expected 
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to be small. As has been said, these predictions could be tested in acceler- 
ator experiments depending on whether the decaying right-handed neutrino 
can be identified. We have discussed how this could be achieved by either 
counting the numbers of leptons emerging from a given vertex or by looking 
to the kinematic endpoint of the lepton pair invariant mass distribution [ID] . 

We have found that the lightest sterile neutrino is a WDM particle which, 
though stable, can not be the only DM component of the Universe. Its 
contribution to the DM relic density is less than 10%. Therefore, other cold 
DM relic must be responsible of the remaining 90%. 

Finally we have derived upper bounds on the Yukawa couplings of the 
model from the experimental upper limit on Br{fi 67). We have shown 
that under the assumption of non-hierarchical Yukawa couplings Br{fi — s> 67) 
is found to be larger than 10~^^ for h > 10~^, i.e within the range of near 
future experiments [33] . 
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